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Narrow Line Seyfert 1 galaxies represent a distinct category in the larger family of Type 1 Ac- 
tive Galactic Nuclei. The occurrence of emission line profiles not broader than 2000kms^', 
combined with luminosity estimates comparable to those of their broad line emitting analogues, 
suggests that some extreme property is characterizing the fueling of their central engines. Two 
hypotheses are more commonly considered: on the one hand, it was suggested that the peculiar 
features of NLS 1 s might be interpreted as the consequence of a nearly pole-on view over a ge- 
ometrically flattened source; on the other, there are arguments favoring the interpretation of a 
relatively low mass black hole, accreting at very high rates. Both explanations provide advan- 
tages and drawbacks, but, in spite of the different assumed physics, they agree in identifying 
NLS Is as the crucial test ground for our understanding of the Broad Line Region structure and 
dynamics in AGN. Here we report the results obtained investigating asynchronous optical and 
X-ray spectroscopic observations, respectively extracted from the Sloan Digital Sky Survey and 
the XMM-Newton science archives, of a sample of both broad and narrow line emitting Seyfert 
1 galaxies. Exploiting data collected by the different instruments carried on board XMM, we try 
to identify the various components which sum up in the observed broad band spectra of the two 
classes. We discuss the relation that thermal, non-thermal, line emission and broad band absorp- 
tion components of the X-ray source show with the optical emission line profiles, as interpreted 
in the framework of a composite Broad Emission Line Region structural model. 
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1. Introduction 

After more than 25 years of extensive investigations, Narrow Line Seyfert 1 (NLSl) galaxies 
still pose some of the most intriguing challenges to our understanding of Active Galactic Nuclei 
(AGN). Following their first detections and the subsequent identification as a peculiar Seyfert class 
[12, 28], it was quickly realized that the physics of nuclear activity is somehow carrying out its 
extreme configurations in these objects. Although they were first noted to have the characteristic 
features of Type 1 Seyfert spectra, but with recombination lines only slightly broader than forbidden 
ones, subsequent studies pointed out many peculiar properties that extend well beyond a mere line 
width based distinction. 

Applying the correlation analysis of Boroson & Green [6] among the emission line profiles 
and intensities, NLSl galaxies were found to posses the smallest values of the [O III] / Hj8 ratio 
and the strongest Fe II multiplets observed in Seyfert Is [33, 34, 7]. Significant line shifts and 
asymmetries were detected in C IV [36, 23, 35] and [O III] [38, 4, 8], while a systematically high 
[S II] 6716/6731 intensity ratio indicates a low density Narrow Line Region (NLR) environment 
[37]. At X-ray energies, the spectra are usually characterized by a steep power-law continuum in 
the hard band (2 - 10 keV) and by an enhanced variability with respect to broad line Seyfert Is (e. g. 
[22, 15, 24]). The optical variability, on the contrary, is not exceptional [13, 16, 3]. Komossa et al. 
[17] pointed out that NLSls are rarely radio-loud, though signatures of jet activity have sometimes 
been spotted [14, 1,2]. 

Several interpretations have been proposed, in order to explain the characteristic features men- 
tioned above. Geometric effects, such as flattening and inclination, as well as winds and density 
properties, metallicity, absorptions, and high accretion rates onto relatively low mass black holes 
were taken into account (see [18] for a comprehensive review). At present, the most widely ac- 
cepted paradigm assumes that NLS 1 properties mainly descend from the small mass of the black 
hole powering the central engine, estimated to be in the range of 10^ — IO^Mq, against the typical 
masses of 10*^ - 1O'^M0 found in broad line Seyfert Is (BLSls) and QSOs. Since both NLSl and 
BLSl nuclei span the same range in luminosity [25], it follows that the accretion rates of NLSls 
are considerably higher than in BLSls. This interpretation agrees with the observed relationship 
among the black hole mass and the luminosity of the host galaxy bulge, which is generally faint in 
NLS Is [9], but it raised some problems in the corresponding relationship associated with the stellar 
velocity dispersion. Due to the difficulties in getting accurate determinations of stellar kinematics 
in AGNs, however, the velocity dispersion of the narrow line region gas has been sometimes used 
as a surrogate. Botte et al. [10] questioned the validity of such a proxy, pointing out that additional 
contributions to the gas kinematics, driven by winds associated to the high accretion rates of the 
central engine, could originate the discrepancy, which, instead, is appreciably reduced when direct 
measurements of the stellar velocity dispersion are available. 

As a contribution to the ongoing investigation on the nature of NLSl galaxies, here we present 
a compared analysis of the optical and X-ray spectra of BLS 1 s and NLS Is. We study the role played 
by geometric effects and intrinsic properties of the central engine in the formation of the observed 
emission line profiles and their influence on the ionization conditions of the line emitting plasma. 
For this reason we shall organize our discussion as follows: in §2 we present the relationship 
among broad line plasma physics and the observed emission line profiles; in §3 we describe the 
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observations taken into account in our analysis, while §4 provides an overview of our results and 
the final discussion. 

2. Plasma diagnostics with Boltzmann Plots 

Type 1 AGN spectra are fundamentally characterized by the presence of a strong non-stellar 
continuum and by several broad and narrow emission lines. While it is observed that a narrow 
component appears in every emission Une, the broad components, on the contrary, mostly affect 
the permitted and some inter-combination lines. It is, therefore, clear that the origin of line emis- 
sion can be traced back to different regions, having distinct physical properties. In the so called 
Narrow Line Region (NLR), which is located at a large distance from the central continuum source, 
so that the emission lines have profiles with FWHM < 10^kms^\ we find a low density plasma 
(Ne ~ 10^ cm^^), where the meta-stable levels of ion species, which are charged by collisional 
processes, can survive long enough to fall back to their ground configuration with spontaneous 
radiative decays. Deep in the gravitational potential of the central engine, much closer to the con- 
tinuum source (at a distance R < I pc), instead, the Broad Line Region (BLR) hosts a high density 
medium (10^ cm^^ <Ne< lO^^cm"^), where forbidden lines are suppressed by the collisional 
de-excitation of their upper levels. Only the emission fines with high transition probabilities can 
originate in the BLR and the kinematics of the line emitting gas results in a profile broadening with 
FWHM > lO^kms"^ 

The suppression of forbidden lines, as a consequence of the significant influence of collisional 
processes, makes it difficult to estimate the BLR physical properties, by means of the usual diag- 
nostic techniques applied to other nebular environments. However, if we consider an emission line, 
originated by the transition from an upper level m to a lower level Z, we can express the associated 
intensity as: 



where A„/ is the probability of a spontaneous radiative decay, Nu{r) the number density of the 
emitting ions, x{r) is the optical depth of the layer, h, c, and A„/ represent the Planck constant, the 
speed of light, and the line wavelength, while R\ and R2 are the boundaries of the emitting region. 
In the case of an optically thin plasma, with the population of the high excitation energy levels 
following the Saha-Boltzmann distribution (plasma in a Partial Local Thermodynamic Equilibrium, 
or PLTE; see [26, 27]), Eq. (1) takes the form of: 



in which gu is the statistical weight of the upper level, t is the spatial extension of the emitting 
region. No is the number density of the ions in their unexcited configuration, while £",,0, kp, and T 
are the excitation energy of the level, the Boltzmann constant, and the plasma electron temperature, 
respectively. Introducing a normalized line intensity, with respect to the atomic constants that 
characterize the transitions, we find from Eq. (2): 




(1) 
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Figure 1: Four examples of Boltzmann Plot applied to the broad Balmer emission line components: (upper 
left) no linear fit to the Balmer series; (upper right) poor fit to the expected linear behaviour; (lower left) 
the BP provides an appreciable fit to the Balmer lines up to H5; (lower right) the broad components of the 
Balmer series are detected up to He and they fit in the linear function. 

If we take into account a series of transitions originated by a particular ion species, the normalized 
intensity of each hne is a function of the upper level energy, that we can express as: 

logg 

log/,, = --—EuQi + const., (4) 
Kb! 

where it has been assumed that the line emitting region is the same for the whole transition series. 

Eq. (4) defines the Boltzmann Plot (BP) of the transition series, which, in the framework of 
the assumptions described above, predicts the normalized line intensities to behave as a specific 
function of their excitation energies and of the plasma electron temperature. The BP method, 
therefore, represents a valuable tool to investigate the physical conditions in a line emitting plasma, 
where coUisional processes are frequent enough to affect the population of excited ion species. In 
La Mura et al. [20], this idea was explored with the application of the technique to a statistically 
relevant sample of type 1 AGNs, extracted from the data release of the Sloan Digital Sky Survey 
(SDSS). As it is illustrated in Fig. 1, where some examples of BP applied to the broad components 
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Figure 2: The temperature parameter A plotted vs. the FWHM (upper panel) and FWZI (lower panel) of 
broad emission line components for the different BP results: asterisks are objects with no reasonable fit to 
Eq. (4); filled triangles are objects with a poor linear fit; open squares represent sources with a good fit and 
the broad emission line components detected up to H5; open circles are objects featuring a good fit and the 
broad emission fine component up to He. Data from [20]. 



of the Balmer lines are given, the technique led to the identification of four main possibilities, 
corresponding to different physical scenarios in the line emitting plasma. 

In general, it is observed that a satisfactory match with the underlying assumptions is achieved 
in approximately 30% of the selected objects, with an improved statistics in the domain of sources 
with very broad emission lines. This result is illustrated in Fig. 2, where the temperature parameter 
A = —loge/kBT, corresponding to the slope of the linear function defined in Eq. (4), is compared 
to the width of the broad emission line components. In the range of narrow line emitting objects, 
instead, the BP finds indications of a higher degree of plasma ionization, probably due to a stronger 
interaction with the ionizing radiation field, which leads to the gradual loss of the equilibrium con- 
ditions required by the method. Furthermore, in the objects where the fundamental BP assumptions 
are still appreciably holding, the inferred temperatures are averagely higher than the values found 
in broad line emitting sources. 
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Table 1: List of targets selected from the combined optical and X-ray observation sample providing the 
object names, coordinates (J2000), measured X-ray fluxes, and redshifts. 
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IV. r\. 




rx JVC V - iz Kc V J 
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(h:m:s) 


(deg:m:s) 
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0.10/ 
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09:23:43.01 


-1-22:54:32.4 
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11:1 /:06.39 


-1-44:13:33.3 


o c n 1 _L A on 
35.91 ± 0.29 


A 1 /I /I 

0. 144 


Ti/~^ 1 1 1 c 1 /I f\n 


1 1 .1 CIJA OA 

11:18:30.30 


1 /IA.OC.C/1 1 

-1-40:25:54.1 


4l.lv ± 0.34 


A 1 CC 

0.155 


zjj/iZiD.y+u /uu 


iz. ly. ju.o / 




1 Q ^; 1 4- A AO 
lo.Dl ± U.+Z 


U.Uol 


Mrk50 


12:23:24.14 


-^02:40:44.4 


278.17 ±4.70 


0.023 


Was 61 


12:42:10.60 


-^33:17:03.0 


131.11 ±0.29 


0.044 


PG 1352-H183 


13:54:35.69 


-^18:05:17.5 


55.99 ±0.55 


0.151 


Mrk 464 


13:55:53.52 


+38:34:28.7 


47.17 ±0.75 


0.051 


PG 1415-^451 


14:17:00.83 


+44:56:06.3 


33.49 ± 0.26 


0.114 


NGC 5683 


14:34:52.48 


+48:39:42.9 


118.25 ± 1.50 


0.036 


Mrk 290 


15:35:52.42 


+57:54:09.5 


173.49 ±0.71 


0.030 


Mrk 493 


15:59:09.67 


+35:01:47.3 


140.22 ± 0.69 


0.031 



3. Data description and analysis 

The detection of larger temperatures or, more generally, of an increasing degree of ionization 
in the BLR plasma of objects with narrow emission line profiles is an intriguing point, suggest- 
ing that the characteristics of the ionizing radiation observed in NLS 1 galaxies are different with 
respect to their broad line emitting counterparts. Such differences can be naturally expected as a 
direct consequence of the high accretion rates required by the interpretation of NLS Is as low mass 
objects. However, the gradual loss of fundamental BP assumptions in the range of strongly pho- 
toionized plasmas forces some care in the interpretation of the results. An independent test of the 
observations involves the study of type 1 AGN spectra over a broader range of frequencies. 

Here we report on the analysis of a sample of objects, for which high quality spectroscopic 
information could be retrieved by matching the SDSS Data Release 7 database with the XMM 
Newton observation archive. A selection of type 1 AGN optical spectra, from objects corresponding 
to X-ray sources brighter than 5 • lO-'^ergcm-^s-^ in the 0.2 keV - 12 keV energy range at the 
time of observation (a flux limit based on the quality of observations with the XMM instruments in 
exposure times of a few 10^ s) led to the identification of 122 potential targets, having both broad 
and narrow emission lines. Work on this sample is still in progress, but Table 1 provides a list of 
the sources that have been studied so far. In the following we describe the data analysis and the 
information extracted from the optical and X-ray spectra. 

3.1 Optical spectra 

The identification of the BLR component in the spectra of type 1 AGN requires a careful data 
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Figure 3: Multiple gaussian decomposition of Hj3 and [O III] XX 4959,5007 from the spectrum of 
MCG+04-22-042. The subtraction of the NLR contributions is performed assuming the well defined [O III] 
X 5007 as a template for the profiles of narrow emission lines. 



handling process, which must account for the several contributions that overlap on the nuclear 
source in the observation. A great advantage of the SDSS database is that it provides a standard 
calibration of data in physical units, therefore largely simplifying the reduction process. The fixed 
fiber spectrograph aperture, on the other hand, implies that various fractions of the host galaxy light 
are recorded, together with the AGN signal, depending on the target redshift. Moreover, extinction 
from the inter stellar medium (ISM), arising in our own Galaxy, must also be taken into account. 

The preliminary reduction process, producing the spectra of our selected targets for subsequent 
analysis, involved the Galactic extinction correction and the subtraction of the signals arising from 
outside the BLR itself. At first, we applied the correction for Galactic extinction, using a continuous 
extinction curve, in the form of the function proposed by Cardelli et al. [11], deriving the absorption 
coefficient Ay in the direction of observation [31]. After the corrected spectra were taken to their 
rest frame, using the peak of the [O III] A 5007 as a reference, we proceeded with the removal of the 
underlying continuum and of the emission line blends, affecting the profiles of the broad spectral 
features. An example of this final step is shown in Fig. 3, where the profiles of the emission lines 
in the spectral region of H/3 are modeled by means of multiple gaussian fits. In this technique, 
we assume that the narrow emission lines can be represented by a standard template profile, built 
to match the [O III] A 5007 spectral feature and scaled to the appropriate intensity, and we use 
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multiple gaussian components to fit the broad contributions. The result is a fair subtraction of the 
NLR spectrum from the observation, but the isolated BLR component still retains some relevant 
concerns. Indeed, it has been pointed out that the broad emission lines are often blended together, 
affected by the contamination of the Fe II multiplets, and probably complex in their own origin. 

In this report we shall not deal with the complete treatment of the techniques that are required 
to properly account for the complexity of the broad line profiles and for the problem related to the 
Fe II multiplets, but we note that, while the former is usually a matter of debate in the domain 
of broad line emitting sources, the latter is of particular importance for NLSls. Here we limit 
our analysis to the identification of those contributions that are most likely arising as an emission 
within the BLR and to the measurement of their fundamental properties. The interested reader may 
find more details on these topics in the literature (see e. g. [21, 19]). 

3.2 X-ray spectra 

In the X-ray domain, the energy resolution of our instruments is in most cases still too low for 
a proper detection of spectral line features. High dispersion X-ray spectroscopy in the low energy 
band (0. 1 keV <E< 2.0keV) is possible in the case of bright objects, but it requires large amounts 
of observational time to provide significant data. Broad band spectroscopy in the energy range 
0.2keV <E< 12.0keV, on the other hand, is much more affordable, although the interpretation of 
results relies onto a different approach, with respect to the classical spectral analysis. Basically it 
is assumed that the observed spectrum can be regarded as the combination of a set of fundamental 
components. These build a global model, which, after the instrument characteristics have been 
taken into account, is used to fit the observation and to estimate the identification likelihood of 
every component. 

The X-ray data considered in our work have been collected with the MOSl, M0S2, and pn 
detectors of the EPIC camera caiTied by the XMM observatory. The calibration of observations 
was accomplished by means of the standard procedures of the XMM Science Analysis Software 
(SAS version lO.O.O), using the emproc and the epproc tasks to extract observations from the 
data packages, then the evselect and tabgtigen to filter the observations from spurious sig- 
nals and background effects. We collected our source and background spectra in circular regions 
of 30" in radius and we used the rmf gen and arf gen tasks to derive the instrument response 
information. Finally, the spectra were binned to have at least 15 counts per channel. 

Our choice to use the different detectors of the EPIC camera results in the extraction of approx- 
imately simultaneous observations of the targets, performed by instruments with different energy 
responses. This circumstance, when available, took us to exploit the technique of simultaneous 
modeling of the spectra, where it is required that a single model fits distinct data sets. We give an 
example of this technique in Fig. 4. To build the models, we used the XSPEC software (version 
I2.6.0q), fitting the spectra in several subsequent steps. At first, we ignored the low energy band 
data and we reproduced the spectrum between 2 keV and 10 keV by means of a simple power 
law, combined to a gaussian function at a reference energy of 6.4 keV, accounting for the Fe Ka 
emission line. As a second step, we froze the power law component and we included in the fit the 
low energy band between 0.3 keV and 2 keV. This is a paiticulai^ly complex region of the spectrum, 
where the contribution of unresolved emission and absorption features from highly ionized species 
can be relevant [5]. Since the low spectral resolution and the noise fluctuations did not give us 
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Figure 4: X-ray spectra of MCG+04-22-042 recorded by the MOSl (black), MOS2 (red), and pn (green) 
detectors of the EPIC camera. The continuous Unes represent the model that is assumed to simultaneously 
fit the observations, producing the residual pattern illustrated in the bottom panel. 

the opportunity to significantly distinguisli among tlie possible theoretical models that have been 
developed to account for the spectral energy distribution in this band, we added a combination of 
thermal components, using up to two contributions, depending on the resulting improvement of the 
fit. Once every component had a steady role in subsequent iterations, we allowed for the fit of the 
whole spectrum. The final model was required to account for the X-ray absorption expected from 
the neutral gas column density along the observation direction within our Galaxy. 

In order to estimate the reliability of every component, we only accepted additional contribu- 
tions leading to significant improvements of the fit quality, for which we required a reduced ^ 1 
and a probability of getting such residuals from a correct model not below 30%. Uncertainty es- 
timates on the model parameters were computed at a 95% confidence range by applying the error 
analysis based on the exploration of the parameter space about the best fit model. 

4. Results and discussion 

The purpose of our analysis is to combine optical and X-ray spectroscopic observations in 
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Table 2: Spectral properties measured in the optical and X-ray domains of the considered sample. For every 
object, we provide the FWHM of H/3, the central energy Epe of the Fe emission, the equivalent width EWpe, 
and the power law photon index F, estimated in the 2 keV - 10 keV energy band. 



Name 


FWHMh/3 




EWpe 


r 




(kms-i) 


(keV) 


(eV) 




2MASX J0306395 8+0003426 


1949 ± 170 


6.59 ±0.18 


202.43 ± 88.31 


1.87 ± 0.02 


MCG+04-22-042 


1946 ± 170 


6.52 ± 0.06 


95.55 ±48.14 


1.74 ± 0.04 


Mrk 110 


2339 ± 348 


6.45 ± 0.04 


59.57 ± 16.68 


1.71 ±0.01 


PG 1114+445 


5090 ± 462 


6.42 ± 0.05 


115.47 ±31.05 


1.53 ± 0.05 


PG 1115+407 


1760 ± 170 


7.18 ±0.10 


212.75 ± 130.91 


2.23 ± 0.06 


2E12 16+0700 


2028 ± 170 


6.48 ± 0.42 


412.80 ± 276.98 


2.21 ± 0.05 


Mrk 50 


5540 ± 199 


6.37 ± 0.06 


61.53 ±40.96 


1.79 ± 0.03 


Was 61 


1591 ± 170 


6.56 ±0.11 


95.86 ± 42.68 


2.03 ± 0.03 


PG 1352+183 


4458 ± 829 


6.45 ± 0.36 


179.48 ± 116.27 


1.92 ± 0.06 


Mrk 464 


6188 ± 575 


6.43 ±0.18 


192.70 ± 117.88 


1.54 ± 0.06 


PG 1415+451 


2850 ± 170 


6.67 ± 0.07 


163.30 ±92.51 


1.97 ± 0.05 


NGC 5683 


5273 ±313 


6.41 ± 0.20 


514.44 ± 336.87 


1.74 ±0.15 


Mrk 290 


4202 ± 170 


6.37 ± 0.06 


74.51 ±35.30 


1.54 ± 0.02 


Mrk 493 


988 ± 170 


6.86 ±0.19 


239.25 ± 155.68 


2.10 ±0.05 



order to explore the physics in the BLR environment of AGNs. In Tab. 2 we collect our estimates 
of some fundamental parameters, which we extracted from the optical and X-ray spectra. The most 
notable effects, which we also illustrate in Fig. 5, involve the relationships of the Hj3 FWHM with 
the Fe emission line central energies and the hard X-ray power law photon index. As it was already 
noted [32], we detect a clear increase in the hard X-ray power law slope moving from broad to 
narrow line emitters. On the other hand, reflection of the hard X-ray continuum very close to the 
central source originates a blend of Fe lines, with a resulting energy that depends on the ionization 
degree of the medium [30] . The most prominent contribution to this component descends from the 
Fe Ka emission of the varously ionized Fe atoms. Its energy shows a steep raise in the domain of 
NLSls, indicating a dominant role of ionized reflectors in these objects. Actually, indications of 
reflection within a significantly ionized medium in NLS Is have been previously found (see [29]). In 
this work, however, we used the Fe line central energy to investigate the properties of the reflecting 
medium within a sample of targets including both narrow and broad line emitters. The results that 
we achieved so far suggest that, while in the range of very broad line emitting sources we observe 
radiation coming from a neutral component of the medium, the degree of ionization averagely 
increases when moving towards objects characterized by narrower emission line profiles, with a 
sharp raise in NLSls. 

The detection of an averagely higher ionization degree among NLSls, with respect to broad 
line emitting objects, together with the different properties of the hard X-ray continuum, supports 
the conclusions anticipated by the BP studies of large samples of type 1 AGN optical spectra. In- 
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2000 4000 6000 2000 4000 6000 



FWHMh^ (km s-') FWHMjj^ (km s"') 

Figure 5: Hard X-ray power law continuum photon index F (left panel) and the Fe line central energy (right 
panel) plotted as functions of FWHM[jj3 . Notice the smooth increase of the power law slope, while moving 
from the domain of BLS 1 s towards NLS 1 s, as opposed to the sharp raise observed in the Fe emission energy, 
suggesting that the BLS Is emission involves reflection from a neutral medium. 



deed, the high temperatures or, more generally, the stronger effects of the ionizing radiation field on 
the BLR plasma detected by the method are in appreciable agreement with the results observed in 
the X-ray domain, thus reinforcing the highly accreting low mass black hole hypothesis supported 
by that technique. We should note, however, that none of the properties under examination here 
shows a clear break distinguishing the targets simply on the basis of the width of their emission line 
profiles. Indeed, the complex balance of processes, resulting in the formation of the observed op- 
tical line profiles, may well affect our interpretation of the intrinsic source properties, introducing 
some degree of uncertainty. In La Mura et al. [21], the question of what kind of uncertainties might 
be expected in the determination of the central engine physical properties by means of emission 
line profile analysis was addressed exploiting a complex BLR kinematical model, developed to 
account for geometric effects, such as source flattening and inclination. By reproducing the whole 
shape of the emission line profiles, it was concluded that, although geometry can be expected to 
play a significant role in the formation of the observed spectra, it is not a critical issue for NLS Is, 
with respect to the remaining objects. 

From the results of this work, we cannot yet draw any conclusive remark concerning many 
aspects of the nature of NLS Is, but we notice that the use of multiple frequency observations led 
us to identify consistent results from very different investigation techniques. The contribution of 
this analysis, as soon as it will be possible to study the whole available sample and to perform fur- 
ther observations, however, shall certainly lead to a better view of the intrinsic properties of these 
intriguing objects and to a more detailed interpretation of the innermost structure of type 1 AGNs 
in general. 
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